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We calculate analytically the phase diagram of a two-dimensional square crystal and its wrapped 
version with defects under external homogeneous stress as a function of temperature using a simple 
elastic lattice model that allows for defect formation. The temperature dependence turns out to 
be very weak. The results are relevant for recent stress experiments on carbon nanotubes. Under 
increasing stress, we find a crossover regime which we identify with a cracking transition that is 
almost independent of temperature. Furthermore, we find an almost stress-independent melting 
point. In addition, we derive an enhanced ductility with relative strains before cracking between 
200-400%, in agreement with carbon nanotube experiments. The specific values depend on the 
Poisson ratio and the angle between the external force and the crystal axes. We give arguments 
that the results for carbon nanotubes are not much different to the wrapped square crystal. 

PACS numbers: 62.20.F-, 61.46.Fg, 61.72.Lk, 64.70.dm 



I. INTRODUCTION 

The discovery of macroscopical two-dimensional (2D) 
graphene sheets by mechanical cleavaging [J] has demon- 
strated that free-standing or suspended 2D crystals can 
exist despite their large 2D positional fluctuations. Since 
then, a variety of other free-standing 2D crystallites have 
been prepared 0. These crystals are stabilized by fluc- 
tuations in the solid plane as verified experimentally in 
Ref. d, following the predictions in Ref. |3- The wrapped 
version of the 2D free-standing graphene had been found 
much earlier in 1991 Due to their high strength, the 
mechanical properties of such materials have recently at- 
tracted great interest. 

The behavior of three-dimensional (3D) crystals as 
a function of stress is well known. For small stresses, 
they expand elastically with a linear stress-strain curve. 
Above the yield point the curve flattens due to the irre- 
versible plastic deformation. At even higher stress, cleav- 
age sets in with further fracture. If the plastic region is 
small or absent, the material is called brittle, otherwise 
ductile. 

A similar stress-strain curve was expected for 2D crys- 
tals or their wrapped versions. Yakobson et al. @ was 
one of the first to determine the cracking strain of single- 
wall carbon nanotubes (SWNT) by computer simulation. 
Since then, there have been many similar studies us- 
ing different simulation methods (see 0] and references 
therein). Most of these observed a cracking strain be- 
tween 15% and 40% depending on the chirality of the 
tube at room temperature T « 300K. The results of the 
different simulations differ widely. Experimental values 
for ropes of SWNTs found cracking strains of 6% [Sl, l9 | 
and 13% [13] for multiwall nanotubes. Huang et al. |ll| 
measured less than 15% for tensile failure at 300K which 



is defined either by the yielding strain for ductile nan- 
otubes or by the cracking strain for brittle ones. At high 
temperatures they were able to go to extreme elonga- 
tions of 270% before cracking. Due to the large temper- 
ature, the SWNTs show an extremely ductile behavior 
with kink motion along the tubes. These were inter- 
preted as defects which do not only perform glide but also 
climb motion in the SWNT [13, [IS] at high-temperatures. 
In Ref. [3, a molecular dynamic simulation for various 
SWNTs at high temperatures was carried out exhibiting 
large defect formations before cracking. 

Due to large activation barriers, the strain value of the 
yield point is mostly dominated by the creation of closely 
lying defect pairs of opposite " charge" forming dipoles at 
low temperatures. The poles separate by glide motions 
at increasing temperatures. By calculating the energies 
of Stone- Wales defects (SW) one finds, by simulation a 
strain value between 6% (arm-chair tube) and 12% (zig- 
zag) [III, [H, [13, where these defects possess negative 
formation energy. The results were consistent with ex- 
periments based on measuring electronic scattering in the 
tube [l^l ■ Plastic behavior for various SWNTs is seen to 
set in at relative strains between 5% and 10%. 

From numerical simulations, we know that the activa- 
tion barriers for S W defects are quite large [l3, HO, HH . 
It depends on the time duration of stress or heat, how 
the SW defects form. For example, brittle SWNTs show 
a defect formation that leads immediately to cleavaging, 
and subsequent cracking. 

The purpose of this note is to study these processes 
with the help of an extension of a model introduced in 
Refs. [2^ . [23 1 to describe crystal melting of 2D and 3D 
lattices. The model contains linear elastic forces coupled 
minimally to an additional integer- valued plastic field to 
allow for defect formation. This model is here extended 
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by an external stress. By modifying this we can inves- 
tigate phase transitions and instabihties of a 2D crystal 
under stress at finite temperatures. We shall restrict our- 
selves mainly to a square-lattice model and its wrapped 
version, for simplicity. The more realistic triangular and 
honeycomb lattices will be treated in the future by ex- 
tending the corresponding melting models [2^. The ba- 
sic physics will not be much different for different lattice 
symmetries. We come back to this point in Sect. VI below 
where we give more arguments that also quantitatively 
it should not change much going from square crystals to 
triangular and honeycomb lattices and its wrapped ver- 
sions. 

The main advantage of square lattices is that one can 
easily calculate partition functions, which require sums 
over integer- valued defect fields. The sums are simplest 
for square crystals with small Poisson ratios [2^ . It turns 
out that SWNTs and graphene are systems with small 
Poisson ratio ~ 0.14, making them well suited for apply- 
ing this technique. The sums can be performed with a 
technique developed for XY models of superfluidity, us- 
ing an so-called inverse Villain (iV) approximation [2^ . 
The defect model in this approximation will briefly be re- 
ferred to as cosine model. In the cosine model, the defect 
aspects can treated by mean-field methods. 

In the following, we will first discuss the phase diagram 
of extended 2D square crystals starting with the phase di- 
agram of the cosine-model in mean-field approximation. 
We shall find a second-order phase transition line which is 
identified as the cracking transition connecting the melt- 
ing point with a point at zero temperature. In addition, 
we encounter a vertical second-order transition line at 
constant temperature starting at the melting point. 

Next, we discuss the full theory without the iV- 
approximation. Wc shall find a similar phase diagram 
where now the second-order cracking line in the iV- 
approximation is almost everywhere a crossover. For 
temperatures near the melting transition our theory give 
relative strains of 200-400% before cracking. This is in 
accordance to the high-strain values of the experiments 
of Huang et al. [ll|, ^ for SWNTs. We find extended 
defect configurations before cracking consisting of homo- 
geneously distributed defect stripes. 

Finally, we will discuss the physics of large wrapped 
square crystals under stress. We find the same phase di- 
agram, stress-strain function and cracking stress for the 
wrapped version of a square crystal as for the 2D ex- 
tended crystal. The main difference lies in the fact that 
for achiral tubes which we define by the property that 
the vector along the circumference of the square tube 
lies not in the direction of a crystal axis show spiral-like 
defect configurations under stress. In accordance to the 
experiments |12| defect glide and climbs are relevant in 



tubes. 

We point out that within our theory it is impossible to 
find the correct yield point at small temperatures where 
plasticity sets in, since our model does not really account 
for the true activation barriers [13, 0, uM ■ Since activa- 



tion energies at high-temperatures are no longer relevant 
because defects overcome the barriers by thermal fluctu- 
ations, we expect that the yielding point tends to zero 
stress leading to an extensive dislocation creep seen in 
the experiments of Huang et al. (llj . This is the tem- 
perature regime where our theory gives the correct phase 
diagram. 

The paper is organized as follows. In Section II we 
state the model. Section HI contains the calculation of 
the phase diagram within the iV-approximation in mean- 
field. In Section IV, we discuss the full crystal Villain 
model. Section V contains a discussion of the true phase 
diagram for a square crystal under stress by taking into 
account the discussions in Section III and IV. We also 
discuss in this section the cracking stress and the relative 
strains as a function of the external stress before cracking. 
In Sect. VI, wc discuss the modifications of our results 
when considering wrapped versions of 2D crystals and 
carbon nanotubes 



II. MODEL 

The partition function used here for the square crys- 
tal was proposed in Ref. [2^ . It can be written in the 
canonical form as a functional integral 



where 
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2 'J 4(A + Ai)W 
- 27ri(Ty {ViU.j -f n,y) \, (2) 
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(3) 



Here vp = is the area of the fundamental cell where 
a is the lattice constant. The exponent in Eq. ([1]) con- 
tains the canonical representation of elastic and plastic 
energies, summed over the lattice sites x of a 2D lat- 
tice. The canonically conjugate variables of the distor- 
tion fields V jUi are the stress fields (Jy for i < j with 
the abbreviation (T21 = The stress field a^j ac- 

counts for external forces applied to the boundary of the 
crystal. The parameter [3 is proportional to the inverse 
temperature, /? = a^^/kBT{27T)^. 

The integer- valued fields nij(x) in Eq. ^ are defect 
gauge fields representing the jumps of the displacements 
field Uj (x) over the Volterra surfaces. The lattice deriva- 
tives Vi and their conjugate counterparts denote lat- 
tice differences for a cubic 2D crystal. In Eq. the 
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defect gauge fields is coupled minimally to the dis- 
placements fields Ui. Note that we do not have this min- 
imally coupling in the stress term H^a [ui\ because the 
external force only acts on the surface of the crystal. 
The measure of functional integration in ([1]) is 



J 'D[u.i,crij,nij]: 
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ij nij(x) = -cx) 



°° da. 



(4) 



where N is the number of lattice sites. 

Let us integrate out the stress fields aij in ([1]). This 
leads to the partition function of an elastic Hamiltonian 
with a minimally defect gauge field under stress . We 
use free boundary conditions for the crystal. These are 
taken into account by separating the displacements field 
integration in the partition function over zero momen- 
tum terms Ui{q = 0) and terms with Ui{q ^ 0) In 
the following we first integrate out the zero momentum 
displacement fields Ui{q = 0). One should now take care 
of this integration due to the following fact: A crystal 
which is homogeneously deformed has three independent 
strain directions instead of two which is suggested by the 
counting of the number of displacement fields. One can 
take care of this by integrating over the three indepen- 
dent strain fields Uij = {ViUj + \7jUi)/2 for q = Q 
instead of the displacement fields Ui{q = 0). This leads 
to the partition function 



(5) 



X exp 



Hd[u,, fly-, rij] + H^o [m,] + H^o^ 



with 

H^o [ui] = VF ^ (Jij (VjUj -I- n^j) , 



Hi. 



E 



8(A + m) 



(6) 
(7) 



The displacements fields Ui(x) of non-zero momentum 
are integrated out in the integration measure ([4]) with 
periodic boundary conditions. Because of this, the first 
term in H^q [ui] is actually zero for homogeneous external 
stress. It is only displayed in ([6]) to exhibit the minimal 
couphng nature of the defect fields in ([5]). The Hamilto- 
nian H^o describes the well-known elastic energy of a 2D 
crystal under a constant stress if no defects are present. 
Due to its similarity with the Villain-model of superflu- 
idity [22I, the model ([5]) will be called the Villain model 
of crystals. 



III. SQUARE CRYSTAL IN COSINE MODEL 

In the following, we restrict our attention to an ex- 
ternal homogeneous stress along the x-axis, i.e. a^j = 
a'^SiiSji. We shall calculate the partition function ^ in 
the iV-approximation in mean-field for v = Q. This was 
done in the case of zero external stress in the textbook 

We now describe the procedure when taking into ac- 
count external stresses. First, we integrate out in ^ 



the stress fields Ui 



Then one can sum in the iV- 



approximation over the integer defect fields by restrict- 
ing the displacements fields to the fundamental cell [2^ . 
This leads up to structural factors being a function of a 
parameter (}y related to the parameter (3 by an inverse 
Villain transformation [231 



= - 



(8) 



where is defined by = h{(3) / Iq{(3) and /o,/i 

are modified Bessel functions of the first kind. Over the 
relevant regime treated in this paper, (3y is roughly pro- 
portional to d [2^.[25t. 

The properties of the model H]) can then be calculated 
approximately from the lattice partition function 
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with the cosine Hamiltonian 



H^^ = - ^ cos [27r(ViU2 + V2W1)] 



2 cos 



27r(ViUi + — ffO) 



2 cos [27rV2U2] . (10) 



The displacements fields satisfy periodic boundary con- 
ditions. The Hamiltonian ([T0|) represents two one- 
dimensional (ID) XY-models which are coupled by the 
first term. This coupling term causes the melting transi- 
tion for cr° = [Hi. 

There are two identical ways to derive a mean-field ap- 
proximation from the partition function (jlOp . Either 
one uses the Bogoliubov variation principle with a trial 
Hamiltonian, or one inserts constraint fields leading to 
the variational mean-field free energy /var per atom in 
the lattice [l^] . In the following we discuss the first way. 
As a trial partition function we use 



a/2 

^0 = 11 J ^^'^^P 



a.i cos 27r 



_Wj(x) 



(11) 



-a/2 



In general a2 for an external stress which is not 

zero. In Ref. ^ one uses the same Ansatz for the trial 
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FIG. 1: We show the phase diagram in the (1//3, (j'^/2/^)-plane 
of a 2D crystal under stress a'^ calculated in the mean-field iV- 
approximation by using (|13|l . (|f 5p and (|16p . The intersection 
of the transition lines and the temperature axis denoted by 
l//3m is the melting point for zero external stress. 



partition function where ai = a2 when the external 
stress is zero. By using Peierls inequality we obtain an 
upper bound for the actual free energy per lattice site 
/fl = -kBT\n[Z^)/N H given by /var with 



/var = -keTl log[/o(ai)] +log[Jo(a2) 



(12) 



+ Pyld{aifld{a2f + Pv'^cos, y—j hiaif 
+ I3y21d{a2f ~ aild{ai) - a2-fd(a2)| - ^^^2^ ^ 



The best approximation for /a is given by the minimum 
of /var with respect to ai, 012. We mention that the first 
term in (|10p corresponds to the third term in ([T^ . 

The analysis of the saddle point equations for p2|) is 
straight forward. We summarize in the following the re- 
sults briefly. We find four distinct solutions characterized 
by Q!i = or 7^ for i = 1,2. Only three of these sad- 
dle point solutions arc minima of /var in certain (T, a^) 
regimes corresponding to phase regions. To be more spe- 
cific, one can show that the saddle point with ai ^ 0, 
a2 = is not a minimum of /var when comparing its free 
energy value with the other saddle point values. One 
obtains the following three different (3v regimes TZi as a 
fimction of ct°: 



Til : py> (3y^,ia°) , 

Tl2 : 1/2 </3y<Py^,{a°), 

Us : 1/2 >(3y. (13) 



The function /Sy ^((t'') is defined by the implicit equations 
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V,b 
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4/d(a) 



cos ( — ) = -Id{a) - \ll{a) . 
The values of ai in the regimes TZi are 
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TZi : Id{ai) = 

Id{a2) = 

7?.2 : ai = , 
7^3 : ai = , 
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Idia2) ' 
a2 = . 



(15) 



The free energies in the various regimes are given by 

7^l : /var = -A;BT|log[/o(ai)]+log[/o(a2)]-Y/d(ai) 

-!^T ( ) Idic^ifld{a2) \ «f(o-°)^ 
2 2 (2 + /rf(ai)2) / 2fi ' 

A:BT|log[/o(a2)] - f^^l^^)} - ^^^^ , 



Ti-2 : /var 

'■ /var = — 
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(16) 



In Fig. 1, we show the phase diagram calculated with 
the help of ^ and (HH). We obtain two 

second-order phase transitions between the three regions 
TZi. The fact that the transition between TZi and TZ2 
is of second-order type can be best seen by using the 
stationarity condition for /var in the saddle point. This 
transition corresponds to the cracking transition. The 
transition line intersect the 1 //3-axis at the melting tran- 
sition point l//3m corresponding to 1//3 ~ 2.85. On the 
low-temperature side, the transition line intersects the 
(T°/2/i-axis at cr'^/2ii — 1/3. We note that in the regime 
TZ2 where ai = we obtain from p2p the mean-field vari- 
ational energy of a ID XY- model. The second-order tran- 
sition between TZ2 and TZ'i corresponds then to the phase 
transition of a ID XY-model obtained by the help of the 
mean-field approximation. In an exact treatment of the 
ID XY-model this transition is of course not existent be- 
ing only an artefact of the mean-field approximation [2^ . 
There is in fact an argument that this phase transition for 
the 2D crystal under stresses is existent because for real 
physical systems there exist a melting transition beyond 
the cracking of the crystal. Note that we cannot get the 
cleavaged cracked state exactly within our model since in 
real physical systems the time scale beyond cracking are 
so long that the thermodynamical average is no longer 
fulfilled for a concrete system. Within our formalism, we 
can only describe the physics of the Gibb's state includ- 
ing a thermodynamical average. In the next section we 



5 



0.6 - 



0.5 



3 a 
<1 I 



0.15/ 
0.5, 

I.O/' .' ! 

; 



2.5/ 



1.5, , 

■ 



2^ 
2|U 



FIG. 2: We show the relative elongation rates Auy /a in the di- 
rection of the external force as a function of the dimensionless 
external stress a'^ /2^. The (black) straight solid curve corre- 
sponds to two times the elastic elongation, i.e. 2cr°/2/i. We 
show the elongation rates for various different dimensionless 
temperatures 1//3 shown as numbers located at the intersec- 
tion point of the corresponding curve with the (black) solid 
straight curve. The x-axis value of the intersection point is 
then given by the dimensionless cracking stress. 



shall discuss the full model H]) with ^ were we still find 
a melting transition independent of the external stress 

Next, we calculate the relative strain parallel to the ex- 
ternal force within the mean-field approximation of the 
cosinc-modcl. It is given by Au||/a = 9/vai/9(vFO'*')- 
In Fig. 2 wc show relative strain values for various di- 
mensionless temperatures 1//3 as a function of the exter- 
nal dimensionless stress ct°/2^. The solid (black) curve 
shows two times the elastic part of the total relative 
strain, i.e. 2CT°/2/i. We do not obtain from the mean- 
field approximation of the cosine-model elongation rates 
before breaking of a few hundred percent seen in carbon 
nanotube experiments [ll| but only less than 2/3. The 
reason lies presumably in the mean-field approximation 
which does not take into account the defect degrees of 
freedom correctly. This will be shown in the next section 
by taking into account the stress degrees of freedom in 
the full crystal Villain model exactly when calculating 
the free energy. 



y — 0. We start by integrating 



in Rcf. '22' for the case ct?, 

out in ([5|) the displacements fields Ui and afterwards the 
stress field cr-y . This leads us to the high-temperature 
limit of the partition function Z 



with 



n 



X exp 



E 

x(x)ez 



"TaE^W" ^(x-x')x(x') 



(17) 



(18) 



where w(x) is short for A ^(x) and A(x) is the 2D lattice 
Laplacian. /3 is defined by /3 = /3(1 + A/(2yU + A)) and 
we shall use further below the abbreviation /i = ^(1 -|- 
A/(2/i-|- A)). The lowest-order approximation to the high- 
temperature partition function Z(f is given by [23| 



(2^/3) 



-3N/2 



/i 



4(A + /i) 



N/2 



(19) 



Taking only into account the dominant stress configura- 
tion x(x) = if^x^xo we obtain [2^ 



^stress « exp[2iVe" 



-5/13] 



(20) 



Note that we obtain an agreement between the free 
energy density in the high-temperature phase TZ^ (|16p 
for the cosine-model with the cr°-part of the high- 
temperature free energy density calculated from (I17p 
which is given by 2H'^a/N. 

Next, we calculate the low-temperature expansion. For 
a^j = const, one can skip the first term in © (we left this 
additional vanishing term only in ([6]) to show the mini- 
mal coupling form and thus the gauge degrees of freedom 
of the defect fields). First, we integrate out the displace- 
ment fields Ui and afterwards the stress-fields cTy . Then 
one obtains the partition function of an elastic crystal 
under stress times a defect dependent term denoted by 

Zdcf 



Zj — -^0 e <t"' Zdof 



(21) 



IV. SQUARE CRYSTAL IN VILLAIN MODEL 

Next, we discuss the partition function ^ of the full 
model (IS]) , (l7|) . The melting line of the 2D square crystal 
under the influence of homogeneous stress will be calcu- 
lated similarly to the stress free system j2^ . [23j . This 
was done by intersecting the high and low-temperature 
expansion of the partition function ([5]). 

First, we calculate the simpler case of the high- 
temperature expansion. This was carried out formerly 



with the lowest-order result [22, u2 



Z, 



T-+0 



(271/3) 



N/2 

^ ^ e-^^ (22) 



A + 2^ 



where £ w 1.14. The defect part of the low-tcmpcrature 
partition function is given by 



^dcf 



E E °^^P 



1 



, (23) 
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with 



= 47r2;3^[e,,,V,ejyV,n,.y(x)] 

x,x' 

X t-2(x - x') [ekk'VUwViUk'ui^')] , (24) 

(25) 



H^o [n] = +VF cr° (x)ri,j (x) 



The symbol S denotes the set of gauge-inequivalent de- 
fect configurations Uij on the lattice which omit all gauge- 
equivalent versions n[j = n,;j(x) -|- V^Aj with periodic 
functions A;. The first term in (j23[) is a pure defect in- 
teraction energy term, the second term a stress-defect 
interaction term. The latter results in the well known 
Peach-Koehler force when calculating the force on a de- 
fect configuration due to the external stress [l^. 

In the following, we choose cr^^ = a'^SnSji as in the last 
section, and consider the following defect configuration 



,(x) = ±S.i^iSj^i Y[ (^x 



(26) 



From (|24p we find that this defect configuration has 
an energy dependence i/dcti^-'l/^sr ~ ln(Z). On the 
other hand the stress-defect energy term i/^o [n'] is 
proportional to the length of the defect line I, i.e. 
H^o\ii}]/kBT ~ ±/. This means that for an infinite crys- 
tal one can construct localized defect configurations with 
arbitrarily small energies. This leads us to the conclusion 
that for high-temperatures where activation energies are 
no longer relevant, the yield point where plasticity sets 
in lies almost at zero stress. 

The defects of the type n} make the defect partition 
sum Zdcf in (j23|) diverge, when performing the sum over 
all defect configurations in the free energy, ordered by 
their geometrical size in the well-known cluster expan- 
sion. This is in contrast to the convergent sum if one 
calculates Z^d for a" = ([23]) due to the fact that large 
defect configurations have generally large defect energies 
Hdef [H; [13 and thus small Boltzmann factors. Mathe- 
matically speaking, one cannot interchange the cluster 
expansion sum of (j23[) with the thermodynamic limit 

iV -> oo for crV 0. 

In order to calculate the free energy, we define the fol- 
lowing defect fields 

n 

(27) 

The integer- valued numbers n,; , di are defined below. We 
now Fourier-transform the defect fields riy (x) leading for 
TV ^ oo to 

n,j(k)=^n,j(x)e-*'' (28) 

X 

= (27r)^ 5i,i Sj,i 5{adiki — 27rni)S{ak2 — 27rri,2) 

"1 ,"2 

+ Si,2 <5j,2 S{aki — 2imi)5{ad2k2 — 27r?T,2) 



In the following, we substitute Uij +nij in Eq. p3|) . 

The aim is to choose n^, di in such a way that the terms 
proportional Sxtr^ Tiy vanish in the exponent of (j23p . 
From (|24p we obtain only non-vanishing contributions 
of the Til, n2-sum in (j28p for ni = n2 = 0. In order 
to obtain the correct result we have to be rather careful 
taking the zero momentum limits in (|24p . This should be 
done for a finite lattice system with vertices. In order 
to get the zero momentum limit of pi]) we should first 
take into account that in ([S]) the integration over the zero 
momentum elongations Ui are excluded since this corre- 
sponds to a translation of the solid. By carrying out the 
integration over the stress fields (Ty we obtain the zero 
momentum limit of Hdcfin] 



Hdcf (q ^ 0) 
knT 



^^"I2(x)ni2(x') 



N ^ 2 

X.X 



+ ^j.nii{x)nii{x') + -n,i{x.)njj{-K') (29) 

With the help of this expression, we can rewrite Z^et 
for iV > 1 as 



Zdcl — ^ 



/ cxp 

5 riij G*S 



1 



knT 



{Hdcf [n + n] + H„a [n + n] ) 



J2 







ksT 


■ exp 



X exp 



A 

2^ 



ni a 
di fl 



ni 
di 



111 
d2 



(30) 



where ni/di are determined by the equations 



(2/i + A) 



di 



A 



n2_ 

d2 

(2^ + A)-i + A-i 

02 d\ 



0. 



(31) 



These equations ensure that all terms proportional to 
'^xcfjUij vanish in the exponent of pop . Solving them 
we obtain for nj/dj with j — 1, 2: 



ni a n2 
d\ 2\i ' d2 2fl 

These values simplify the expression ((30|) to 



Zd: 



E 



J2 



knT 



exp 



Hi 



knT 



(32) 



(33) 



Finally we calculate Zdef for (t° = 0. Taking into ac- 
count only the dominant defect configurations nij(x) £ 

as in Refs. [li, [Hi 

Zdci{<y" = 0)~ exp[2cxp(-6.3^) +4exp(-13.7^)] . 

(34) 
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(a) 



a 



FIG. 3: In (a) we show a defect configuration riij = SuSjifiii 
with nil is defined in (I27p for ni = — I and di — 2. Figure 
(b) shows a dislocation pile up in a glide plane leading by 
merging of dislocations to a cleavage being the starting point 
of cracking. 



From the considerations above we conclude that the 
cracking transition is in fact not a true phase transition 
in the Villain model but rather a crossover. The crossover 
temperature can be obtained by the assumption that the 
defect configurations hn or 7122 cover the whole crystal 
area at cracking, meaning that di = 1. A defect config- 
uration which covers only half of the crystal is shown in 
Fig. 3a, where = Si^iSj^ihn with ni = —1 and di = 2. 
The idea behind this crossover temperature comes from 
the fact that most cracking transition models start with 
a pile up of dislocations within a glide where then by 
merging we obtain a cleavage dislocation being the start 
point of cracking [23] . This is shown in Fig. 3(b). 

From ((32|) and the stability criterium [22| < 1 we 
obtain that the lowest stress configuration where we have 
di = I is given for di = 1 and rti = 1 resulting in a 



cracking stress (t° 



2fl 



^ = 1 



(35) 



Comparing this value with the low-temperature crack- 
ing stress in the mean-field approximation of the cosine- 
model given by crj^/2/i ~ 0.33 seen in Fig. 1 we obtain a 
much higher stress here. The relative strain in the crys- 
tal phase is given by Aw||/a = kBTd\n[Z'^-'°]/d{Nv(T°) 
resulting in 



Auii crO 

= 2 — 

a 2jl 



(36) 



Comparing (|36|) with the elastic part of the relative strain 
(T°/2/i we obtain a factor 2 difference. By using (pS]) with 
((35|) we obtain a relative strain at cracking of 200% or 
Au||/a = 2, respectively. 

Let us now address the question whether the defect 
field configuration of Eq. (P7|) is the only defect con- 



figuration leading to the partition function ((33|) . The an- 
swer is negative. From the above derivation we see that 
any fiij arising from the substitution ni 



explained below (|28p . leads to ([55)1 under the condition 
that nii(k) is only non-zero for fc2 = and n22(k;) for 
fci ~ 0. All these fields will be denoted as deject vac- 
uum. Such riij(x) correspond to defect stripes cover- 
ing the whole width of the crystal where the fiaik = 0) 
values are determined from the condition that the term 
proportional to CFijUij vanishes in the partition func- 
tion. From the thermodynamic point of view none of the 
defect configurations ny (x) which fulfill the above con- 
ditions are preferred. On the other hand we used the 
argument di = 1 for determing the cracking crossover 
stress ([35]) which was justified by the defect-merging pic- 
ture in Fig. 3(b). Now suppose that we have additional 
external conditions in the crystal, for example impurities 
or fixed crystal defects generated during crystal growth, 
which lead to the restriction that defects cannot cover 
the whole width of the system. Then the substitution 
riij —>■ TLij + hij cancels the term OijUij in ([23]) only 
partly. This is so since the exact cancellation relies on 
the fact that the only non-vanishing contribution of (I^B)) 
in (|24[) is given by the ni = n2 = term, which has 
zero momentum. Finite-length defect stripes have this 
property only approximately. These also contribute to 
([M)) . Nevertheless, the partition function is still approxi- 
mated by p3p if these residual terms are suppressed with 



respect to the term <Ji. 



From (|24|) we see that 
([33]) is best fulfilled under the stripe length restriction 
for that defect vacuum fiij which has the largest momen- 
tum region in the vicinity of /ci = 0, where nii(k 7^ 0) is 
almost zero. These defect configurations consists of ho- 
mogeneously distributed defect stripes in perpendicular 
direction where the density in this direction is determined 
from ((32|1 . 

Note that the homogeneity can be only be fulfilled ex- 
actly for stresses cro/2/i G 1/N. The average distance dia 
between the stripes is then given by p2p for Ui = 1. The 
cracking condition di = 1 is again given by the fact that 
the homogeneously distributed defect stripes cover the 
whole crystal area where this criterium is justified by the 
defect- merging picture in Fig. 3 (b). This leads immedi- 
ately to the dimensionless cracking stress ([35]) and strain 

Let us finally remark that the requirement of a homo- 
geneously distributed defect stripe configuration as the 
vacuum is also in accordance with the conception that in 
real crystals the external homogeneous stress should be 
relaxed homogeneously across the area by defects. 

We are now prepared to calculate the melting line 
by intersecting the partition functions of the low- 
temperature ()17p and the high-temperature expansion 
(1^ . The result is 

/3 exp[4 exp(-6.3 /3) -f 8 exp(-13.7 /3)] 

X exp[-4exp(-5//3)] « 0.81 (37) 

independent of cr*^. The phase diagram is shown in Fig. 4. 
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FIG. 4: We show the low and high-temperature intersection 
curve of the free energies given by (|37p ((blue) dashed curve). 
The (black) solid curve denotes the cracking transition line 
determined by 



It displays the intersection line ([37|) of low- and high- 
temperature free energies ((blue) dashed curve) and the 
crossover cracking temperature ([35]) plotted as a (black) 
solid curve. We obtain from Fig. 4 that the melting tem- 
perature is given by 1/(3 « 1.42. The melting line is 
given by a first order transition in the case of a square 
lattice [2^ [2^ . This was also found in Ref. [H by using 
computer simulations. 



V. PHASE DIAGRAM OF 2D-CRYSTALS 

So far we have obtained the phase diagram in Fig. 1 by 
applying a mean-field approximation in the cosine-model 
of crystal defect melting, and the phase diagram in Fig. 4 
from the associated Villain-type model. The main differ- 
ence lies in the fact that in the mean-field cosine-model 
the cracking transition between the phases TZi and 7?.2 is 
a second order phase transition but in the Villain-type 
defect model it is only a cross-over. This and the differ- 
ences in the value of the breaking stresses have their 
origin in calculation of the cosine-model partition func- 
tion in a mean-field approximation. It is a well-known 
phenomenon in many physical systems, especially in low- 
dimensions, that quantum- as well as thermodynamical 
fluctuations can destroy a phase transition which appear 
in a mean-field approximation, leaving only a cross-over. 
We have already mentioned an example for this in Sect. 
Ill with the ID XY- model [2^. Summarizing, we expect 
a similar phase diagram as in Fig. 4 for a real 2D crystal 
under stress. The stress-independent melting tempera- 
ture in this figure is in accordance to the fact that in real 
physical systems there exist a melting transition beyond 
cracking. In a triangular lattice, we expect for the melt- 
ing transition line two nearby Kosterlitz-Thouless tran- 
sitions instead of the first-order transition found in the 
square lattice [H, [2§| . 

In addition, our model was shown to possess a 
temperature-independent cracking transition line in 



Fig. 4 which is in agreement with the experimental de- 
termined cracking stress of 3D graphite in Ref. [30|. In 
that experiment, the cracking stress cr° shows an anoma- 
lous temperature behaviour only in a small temperature 
range just before melting where it starts to increase for 
larger temperatures. Such an increase of the cracking 
stress was also seen in other experiments where this in- 
crease starts even for smaller temperatures [3l| . This 
behaviour is not fully understood yet. That the cracking 
stress has a small temperature dependence is expected 
and should be revealed when going beyond the elastic, 
lowest order gradient expansion approximation used here 
when deriving the Villain lattice defect model. These ap- 
proximations are released by using the iV-approximation 
in Sect. Ill to the Villain lattice defect model [l^l leading 
to Hamilton terms beyond elasticity. This leads us to the 
possibility to determine the physics of real crystal models 
which are not restricted to the elastic, lowest order gradi- 
ent expansion with the following conclusion: We expect 
a small temperature dependence of the cracking stress 
separating the phases TZi with 7^2 for realistic models. 
Nevertheless, the melting temperature transition sepa- 
rating 7^2 and 7^3 should not have a stress dependence 
also for more realistic models. Note that terms gener- 
ated in the Hamiltonian beyond elasticity by using the 
iV-approximation are in general not directly connected 
to higher order terms of a certain real existing crystal. 

In the following, we shall generalize our calculation 
of the last section to crystals with stresses irrespective 
of direction. In order to get the cracking stress one 
has to repeat the calculation of Section IV where now 
fiij in ([27)) contains additionally a term of the form 
"-12(54,1 ^j,2(<5a;i,adi2Z + 42,adi2z)- % Carrying out the cal- 
culation we obtain the same defect part of the partition 
function Z^of as in p3p . The cracking transition is again 
determined by the minimal stress where rii/di = 1 or 
'^■i2/c?i2 = 1, respectively. With this condition, we ob- 
tain for the cracking stress 



2/i (1 + v) 



Min ■ 



sin2(^) - ^ 
cos(i?) sin(i?) 



-1 ^ 



(38) 



where we took into account cr°^ — {cos{i}) , sm{-d)) x 
(cos(7?), sin('!?)). Here -d is the angle between one crystal 
axis and the external force. 

We may now calculate the relative strains orthogo- 
nal to the external force by inserting in ([7]) for a^j an 
additional orthogonal auxiliary stress field in order to 
calculate the orthogonal strains by differentiation. This 
leads to relative strain values parallel to the external force 
Alt II /a and orthogonal to it Au±/a of 



Au 



a 

Au± 



2f! 

2/i 



2/i 



(39) 
(40) 
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FIG. 5: The figure shows /2/i 1)38^ as a function of the angle 
i9 between external force and the crystal axis The numbers 
at the curves denote the Poisson ratio i/. 



We show in Fig. 5 the dimensionless cracking stress 
(T{^/2/i ([55]) for certain Poisson ratios i^. Together with 
we obtain parallel relative strains Au||/a « 200 — 
400% where the concrete value depends on the Poisson 
ratio I' and angle i). We note here once more that the 
relative strains (j39l) . ((40|) are only valid for high temper- 
atures or large time scales such that activation barriers 
are no longer relevant in the crystal. 

These large increases in the strain values at high- 
temperatures are in accordance to observations of Huang 
et al [ill . [T3 | for carbon nanotubcs mentioned in the in- 
troduction of this paper. They found generally relative 
elongations which are at least five times higher at crack- 
ing than the tensile failure value at low temperature [s^ . 



VI. WRAPPED SQUARE CRYSTALS AND 
CARBON NANOTUBES 

In the previous sections, wc have examined the be- 
haviour of large square crystals under stress. In the se- 
quel we shall examine the modifications brought about 
by considering wrapped versions of these. They form 
infinitely long thin tubes with a perimeter much larger 
than the lattice constant. Then the curve of the tube is 
irrelevant and the previous crystal model remains appli- 
cable. 

First, we shall consider chiral square tubes. These are 
defined by the property that the vector along the circum- 
ference of the tube lies in the direction of a crystal axis 
leading to periodic boundary conditions in this direction. 
Due to the periodic boundary conditions along the crys- 
tal axis after integrating out the zero momentum strain 
fields described below ([7]) we obtain that the mean-field 
cosine results of Sect. Ill as well as the results for the 
Villain model of Sects. IV and V are also valid. 
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FIG. 6: We show a spiral-like defect configuration for a achiral 
square tube with ni = — 1 and di = 2 by using the defect con- 
figuration (|27p . This defect configuration is consistent with 
the periodic boundary conditions around the circumferent of 
the tube. 



As described in the introduction section an achiral tube 
is defined by the property that the vector along the cir- 
cumference of the tube lies not in the direction of a crys- 
tal axis. Thus the periodic boundary conditions after 
the integration of the homogeneous strain fields are no 
longer in the direction of a crystal axis but in the cir- 
cumferental direction. The mean-field cosine results of 
Sect. Ill are of course still correct in this case. Also the 
cracking stress ([38]) and the strain-stress relations ([39]), 
(|40|) are still valid. But the relevant defect configurations 
discussed below p6p are no longer homogeneously dis- 
tributed stripes but spiral-like defect stripes. One such 
stripe is shown in Fig. 6. By cutting the tube along the 
axis and projecting it on the plane we can use (|27p with 
n-i = ni2 = 1 where we should further take into account 
ni2 defined in the last section. This generalization be- 
comes relevant for external forces not directed along a 
square crystal axis valid for achiral tubes. The distances 
di and di2 between the stripes are determined by the 
periodic boundary conditions of the achiral tube. The 
number of defect stripes for every sort of defects fiy in 
axial direction is governed by the generalized equations 
of ([5T|) used in Section V if we take into account also the 
defect field ni2. This means that we have to substitute in 
the generalized equation of ([3T|) n.i/di and ^12/^12 by the 
defect density along the axis which is given by the number 
of stripes of a special defect type divided by the number 
of faces along the crystal axis. In general this leads to 
the result that the homogeneous distributed stripes do 
not cover the entire tube length. 

In the uncut tube this defect configuration consists of 
long spiral-like defect stripes whose extension in the ax- 
ial direction is as long as possible and consistent with 
([3T|l . The reason lies in the fact that the substitution 
riij — > riij + fiij in (|23|) leads only to a cancellation of 
the external stress fields in the partition function when 
nii(k) is zero for k2 = (and similar requirements for 
'^22, "-12 as is outlined in Sects. IV and V). Only in this 
case the spiral-like defect configurations exactly cancel 
the external stress field in the partition function (|23p as 
described in Sect. IV. Note that this requirement is ex- 
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actly fulfilled only if a spiral-like defect covers the entire 
tube length. If a spiral-like defect stripe is smaller than 
the tube length where the generalized equation of ((5T|) 
applies as described in the last paragraph, the resulting 
stress term in the partition function (|23|) leads only to 
negligible contributions in the free-energy density for in- 
finite length and large perimeter tubes compared to the 
free energy expressions in (j30p . 

That the spiral-like configurations with the longest un- 
broken defect stripes which fulfill ((5T|) are the most rele- 
vant defect configurations is also obviously by the fact 
that this configuration has lowest energy i/def [f^] +Hg-a [h] 
([24ll . (I25|) . where the zero- momentum part ([29ll is in- 
cluded in this expression. 

Due to the generahzed relations corresponding to ((3T|) . 
the length of the spiral-like defect stripes increase for 
increasing stresses, leading to kink pro pag ation observed 
by Huang et al. for carbon nanotubes [ill [T^ . From the 
definition (|27p and its ni2 generalization in Sect. V we 
deduce further that the spiral-like defect length motion of 
stripes fill and 77.22 a re g lide motions, whereas the stripes 
ni2 move by climbs 2^. Both motions were observed in 
the superelongation experiments of carbon nanotubes at 
high temperatures [l3 | . 

All this discussion leads us to the following scenario: 
When increasing the stress applied to the tube, a spiral- 
like defect for every defect type fii and hi2 becomes 
longer by climb or glide, respectively. When one of the 
defect stripes cover the entire tube length, a new spiral- 
like defect of the same type starts to be formed. This goes 
on up to the point when the tube is covered by defects. 
This is when cracking starts. 

Let us end this section by discuss shortly the differ- 
ence between real world carbon nanotubes and wrapped 
square crystals in this theoretical description. A carbon 
nanotube consist of a wrapped honeycomb lattice with 
two atoms per fundamental cell. By taking into account 
that the energy dispersion in the optical sector of the 
lattice displacements is negligible in comparison to the 
acoustical sector [s^ leads to the result that it is suf- 
ficient to consider minimally coupled integer valued de- 
fect fields only in the clastic Hamiltonian of the acoustic 
sector [13. This then leads to a triangular lattice melt- 
ing model since the honeycomb lattice has a triangular 
Bravais lattice. We have considered in Ref. [H a defect 
melting model for the triangular lattice ([T]) for zero ex- 
ternal stress in the simplest way. The strain-stress re- 
lations ([39]) . (|40|) remain of course still valid for carbon 
nanotubes since different lattice symmetries have no in- 
fluence on the long-range behaviour of the lattice like the 
zero momentum strain-stress relations (|39p . (j40p . Only 
the cracking stress (p8|) and the defect vacuum configura- 
tions are changed. As described above, defect configura- 
tions in a square crystal or its wrapped version consists of 
defect stripes which are directed along the crystal axes 
building in general spirals in achiral tubes. The same 
thing is true for the stress release in carbon nanotubes. 
In contrast to a square lattice where every vertex cuts 



two inequivalent lines along the crystal axes, every ver- 
tex cuts in a triangular lattice three inequivalent lines. 
This leads to the fact that the number of inequivalent 
defect stripes in the triangular lattice is thus a factor 3/2 
larger than in the square lattice. By taking into account 
the defect-merging picture in Fig. 3(b) for cracking we 
conclude that in carbon nanotubes the cracking stress 
should be correspondingly larger than in square tubes. 
The stress release now takes place on a larger amount 
of different homogeneously distributed defect stripes. To 
be more specific, we expect in a first rough approxima- 
tion for carbon nanotubes cracking stresses being in the 
average a factor 3/2 larger than in square tubes. Note 
that this factor is similar to the ratio between the melt- 
ing temperature difference factor of square crystals and 
triangular ones [2^ as well between honeycomb lattices 
and square crystals [3^. In the latter case, one has to 
take into account properly the definition of the elongation 
fields in the elastic acoustic Hamiltonian as a function of 
the atomic elongations [131 in order to obtain a temper- 
ature reduction factor at melting in comparison to the 
triangular lattice. 

A more elaborate treatment of the cracking stress of 
triangular or honeycomb lattices needs much more afford 
which is work in progress. 



VII. SUMMARY 

Motivated by recent experiments revealing [ll| that 
carbon nanotubes under external stress show a strong 
ductile behavior at high temperatures with extremely 
large relative elongations before cracking, we have cal- 
culated in this paper the phase diagram, the cracking 
stress and relative elongations of a 2D square crystal 
lattice. The results are hoped to be applicable to car- 
bon nanotubes, although these form honeycomb lattices. 
By starting from a Villain-type lattice defect model we 
have derived in Section II, using the inverse Villain ap- 
proximation, an XY-likc model for crystals. When cal- 
culating within this model the phase diagram in mean- 
field approximation in the (T, cr°) plane we have obtained 
two phase transition lines (see Fig. 1). We have found 
a second-order transition at lower temperatures which 
we identified as the cracking transition line and a verti- 
cal second-order line as the melting transition beyond 
cracking. The dimensionless cracking stress cr^/2/i as 
well as the relative strain rates parallel to the exter- 
nal force Au||/a have upper bounds of cr[^/2/t < 1/3 and 
Au||/a < 2/3 (see Fig. 2). 

Next, we have calculated the phase diagram of the 
full Villain model within low- and high-tcmpcrature ex- 
pansion of the free energy. Here we have found within 
the low-temperature expansion that the cracking tran- 
sition is in fact not a phase transition but a crossover. 
The crossover line is identified by the requirement that 
the ground state defect configuration under stress should 
cover the whole plane of the crystal. Within our model. 
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the cracking stress is independent of temperature. By 
using the intersection criterion of the low and high- 
temperature expansions of the free energy we obtain a 
melting temperature which is independent of the exter- 
nal stress. The whole phase diagram for external forces 
along one crystal axis is shown in Fig. 4. 

We have deduced in Sect. V from the considerations 
above that a crystal under stress should show, in gen- 
eral, a phase diagram as in Fig. 4. The phase transi- 
tion line in the cosine-model observed within the mean- 
field approximation should vanish upon taking fluctua- 
tions into account, converting it into a cross-over line as 
was shown in Sect. IV within the Villain lattice defect 
model. Nevertheless, we expect in accordance with the 
iV-approximation to the Villain model. Nevertheless, we 
expect by going beyond the elastic, lowest order gradi- 
ent expansion approximation used in the Villain model, 
a small temperature dependence of the cracking stress, 
but no stress dependence of the melting transition tem- 
perature. This is motivated by the results for the cosine 
model of Sect. Ill since the iV- approximation used to de- 
rive this model from the Villain model generate Hamilton 
terms beyond the elastic approximation. 

Finally, we have calculated the cracking stress, and the 
relative elongations before cracking. In Fig. 5 we have 
shown the resulting dimensionless cracking stress as a 



function of the angle between the external force and the 
crystal axes for various Poisson ratios. We have found 
dimensionless cracking stresses cr°/2/2 between 100% and 
200% at high temperatures where potential barriers for 
defects are no longer relevant. The full relative strains 
Au||/a in the direction of the external force are twice as 
large as the elastic strain part. The reason lies in the de- 
fect degrees of freedom which then results in full relative 
strain rates AMy/a of 200 — 400% at breaking depend- 
ing on the direction of the external force and the Pois- 
son ratio. The large difference in the breaking stresses 
and the strain rates between the mean-field result of the 
cosine-model and the exact calculation of the Villain lat- 
tice defect model is presumably due to the mean-field 
approximation. 

In Sect. VI we have obtained that the cracking stress 
relation (p8)) as well as the strain-stress relations (|39p . 
(PO)) are also valid for wrapped square crystals. The de- 
fects are spiral-like for achiral tubes where defect glide 
and climbs are relevant in accordance to experiments. 
For honeycomb lattices or carbon nanotubes also the 
stress-strain relations ([M)) . ([3D]) are fulfilled but the 
cracking stress ([551) now modified. We have argu- 
mcnted in a rough approximation that cracking stresses 
in carbon nanotubes should be in the average a factor 
3/2 larger than in square tubes. 
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